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We discuss the production of pion pairs with a large invariant mass in hadronic collisions. We 
present a study based on a perturbative QCD calculation at full next-to-leading order accu- 
racy, implemented in the computer programme DIPHOX. We provide a comparison of various 
predictions with the corresponding observables measured by the E706 fixed target experiment. 



■ Discrepancies between data and next-to-leading order calculation are carefully analysed. We 

classify the theoretical next-to-leading order distributions with respect to their infra-red sensi- 
tivity, and explain which distributions need improvements. Further, we comment on the energy 
scale dependences of non pertubative effects. 



1 Introduction 



Quantum Chromodynamics (QCD) is now accepted as the theory of the strong interaction. With 
the starting run II of the Fermilab Tevatron, and the forthcoming Large Hadron CoUider (LHC) 
at CERN, QCD is entering a new area associating improved precision measurements with more 
accurate calculations. This concerns not only a refined study of the strong interaction itself and, in 
particular, the ability of perturbative QCD to describe high energy phenomena precisely, but also 
the accurate prediction of sensitivities in Higgs boson and expected new physics searches. Both need 
an accurate understanding of respective signals and backgrounds, which are heavily affected by QCD 
radiative corrections. Meanwhile, the proper understanding of the large amount of data collected by 
fixed target experiments remains rather challenging. In all these respects, the production of neutral 
pion pairs with a large invariant mass is an interesting phenomenon. This process has been recently 
measured by the fixed target E706 Experiment [Q at the Fermilab Tevatron. The study of this 
process and the examinations of the various observables which have been measured are instructive 
for several reasons. 

First of all, it is interesting to compare the ability of perturbative QCD to describe single pion 
and pion pair production respectively. In particular, one can identify the observables for which 
non perturbative effects are negligible and the QCD improved parton model based on coUinear 
factorisation can be safely used, and examine otherwise when it fails and why. In this respect, loosely 
speaking, inclusive enough correlation observables with symmetric selection cuts are less sensitive 
to non perturbative effects such as, e.g. a kx kick, than one particle inclusive distributions 
especially at fixed target energies. Dipion observables thus provide a good opportunity to cross- 
check the consistency of the sets of parton-to-pion fragmentation functions Di^^o{z, M'j), such as 
those given in Ref. |^ and the more recent ones in Ref. Q, which are available in the literature[]. 
Both these sets have been extracted from data on hadronic final states in e^e~ annihilation. They 
are given by analytical ansatze whose parameters were fitted in the range 0.1 < z < 0.9 but which are 
tightly constrained^ only in 0.1 < z < 0.7. In particular, the gluon-to-pion fragmentation functions 
are not much constrained because the subprocesses involving gluons in e+e^ annihilation appear 
at higher order in the perturbative expansion. Instead, subprocesses involving outgoing gluons 
contribute already at lowest order in hadronic collisions. Therefore the contributions from gluon 
fragmentation in hadronic collisions are not negligible and large uncertainties on their magnitude 
may affect quantitative predictions. 

In addition, the production of 7r° pairs, together with the production of pairs 77r°, where 7 is itself 
a prompt photon, is the dominant background against which the experiments measuring prompt 
photon pairs have to fight. The expression "prompt photon" means that such a photon does not 
come from the decays of a hadron such as tt", 77, etc. produced at large transverse momentum. 
This background is especially worrisome in collider experiments, such as the Collider Detector at 
Fermilab (CDF) |^ and DO Q experiments at the Fermilab Tevatron, where the pions cannot be all 
reconstructed event by event. In addition, these two processes of 77r'' and tt'^tt^ production represent 
an important fraction of the so-called reducible background to the search of neutral Higgs bosons 
in the channel — > 77 in the intermediate mass range 80 - 140 GeV at the CERN Large Hadron 
Colhder (LHC) Since the jet-jet cross section is about eight orders of magnitude larger 

than the expected Higgs boson signal for a standard model Higgs boson, this reducible background 
is overwhelming before any selection cut is imposed. Stringent isolation cuts are applied in order 
to reduce this background severely and make prompt photons measurable. Nevertheless, it is a 
theoretical challenge to understand the remaining background at a quantitative level. In this respect, 

^Similar consistency checks for the set determined in especially for the gluon fragmentation functions, have 
been presented in Ref. in the case of hadro- and photo-production of single inclusive charged hadrons. 

^To be precise, regarding pions, the authors of M determine fragmentation functions into charged pions then 
deduce the ones into neutral pions relying on isospin symmetry. 
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the understanding of tt'' pair production in the fixed target energy range wiU aUow one to transport 
this knowledge to coUider energies |Q using the factorization property of perturbative QCD. 

In the present article, we present a phenomenological study of tt^ pair production at Next-to-Leading 
Order (NLO) accuracy in perturbative QCD, based on the computer code DIPHOX In section 
H we briefly describe the theoretical content on which the present study is based. We then compare 
DIPHOX-hased theoretical results with data of the E706 experiment in section ^. Section ^ is then 
dedicated to a detailed phenomenological discussion of the various observables measured in ir'^ pair 
production. The aim of this phenomenological study is 1) to test the ability of the theoretical 
framework used to describe various correlation observables that have been measured by the E706 
experiment, 2) to identify and understand the possible discrepancies between theory and data for 
these observables, and 3) to examine how the sensitivity to non perturbative effects, which limit 
the predictive power of the perturbative approach, evolves as the 'hardness' of the probes increases. 
To this purpose we classify the various observables according to their sensitivity to the infrared 
dynamics. At this occasion, the analysis of the transverse distribution of each pion leads us to 
reconsider briefly the px spectrum in the case of inclusive one pion production as well. Section ^ 
gathers our conclusions and perspectives. The goal is to foresee how reliable predictions based on our 
present knowledge can be for dipion production at much higher, collider energies. In a forthcoming 
paper, we will investigate more specifically our ability to predict the reducible background to the 
Higgs boson search at LHC in the channel H 77, relying on the informations obtained in the 
present study. 



2 Theoretical content 

The computer program DIPHOX on which the present study relies is a Monte Carlo code of partonic 
event generator type. This code has been designed to describe the production of pairs of particles in 
hadronic collisions, accounting for all contributing partonic processes at full NLO accuracy. These 
particles can be prompt photons or hadrons, in particular neutral pions. This code is fiexible 
enough to accommodate various kinematic or calorimctric cuts. Especially, it allows to compute 
cross sections for isolated jtt'^ and tt'^tt'^ pairs, for any infrared and coUinear safe isolation criterion 
which can be implemented at the partonic level. The physical content and schematic description of 
DIPHOX has been given in , where we have concentrated on the case of diphoton production^ 
|lj,|l5|. However, the production mechanisms of fragmentation type involved in diphoton production 
are exactly the same as the ones at work in the hadroproduction of 7 7r° and Tr^Tr". Let us remind 
this point briefiy. 

A prompt photon may be produced according to two possible mechanisms: either it takes part 
directly to the hard subprocess (let us call it call "direct"), or it results from the coUinear frag- 
mentation of a parton which is itself produced with a large transverse momentum (let us call it 
"from fragmentation"). Consequently, pairs of prompt photons with a large invariant mass may 
be produced according to three possible mechanisms: both of them take part directly to the hard 
subprocess (direct mechanism) , or one of them takes part directly to the hard subprocess while the 
other one results from the fragmentation of a large transverse momentum parton (one-fragmentation 
mechanism), or else each of them results from fragmentation (two-fragmentation mechanism). For 
a general discussion of the production mechanisms of prompt photon pairs, see [^S) . Similarly, 7 7r° 
pairs with a large invariant mass may be produced according to two possible mechanisms: either 
the prompt photon takes part directly to the hard subprocess (direct mechanism) , or it results from 
the coUinear fragmentation of a parton which is itself produced with a large transverse momentum 
(fragmentation mechanism) . 

^Thc acronym DIPHOX stands for "diphoton cross sections". 
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In agreement with the factorisation property - even in the case of reasonable isolation requirements, 
see Ref. [ p^ - the partonic subprocesses involved by the one-fragmentation mechanism in prompt 
photon pair production are the same as the ones involved in 7 tt^ production when the prompt photon 
is of the 'direct type'. The only change is the fragmentation functions involved: fragmentation into a 
photon in the first case, into a pion in the last one. Similarly, the partonic subprocesses involved by 
the so-called two-fragmentation mechanism in prompt photon pair production are the same as both 
the ones involved in jtt^ production when the prompt photon come from fragmentation mechanism, 
and those involved in tt^ pair production. Therefore, by keeping only the relevant mechanisms, and 
replacing the fragmentation functions appropriately, DIPHOX allows also to study the production 
of 7 tt'^ as well as of tt*' pairs with a large invariant mass and individual transverse momenta, at full 
NLO accuracy. 

3 Comparison with the E706 data 

The E706 experiment |0] studied collisions^ on various fixed targets (H,Be,Cu) of several beams. 
They used proton beams at different beam energies, 530 GeV {VS ~ 31.55 GeV) and 800 GeV 
{VS ~ 38.76 GeV), as well as a 515 GeV tt" beam (y/S ~ 31.10 GeV). We have decided to focus 
on comparisons with data from proton initiated reactions, because the sets of available pion's pdf 
are rather old and not as reliable as the recent ones for protons. We have also focused on data on 
Hydrogen and Beryllium targets, because the Beryllium nucleus is still rather small hence nuclear 
effects are hopefully reduced, in contrast to the case of a Copper target where nuclear effects are 
more important. We have performed a comparison between our NLO QCD results and data on 
Hydrogen and Beryllium targets, for protons beam energies of 530 GeV and 800 GeV. We have 
chosen to present only the comparison with the data obtained in proton-Beryllium for which the 
statistics is larger than for Hydrogen target. We restrict the presentation mainly to a comparison 
with data from the collision of the 530 GeV proton beam, in order to limit the proliferation of plots 
in the article. Similar conclusions have been obtained in the other cases studied. 

The NLO QCD predictions presented here are provided by the computer code DIPHOX presented in 
the previous section. Both the theoretical results and the data are expressed 'per effective nucleon' 
in the Beryllium nucleus {Z=4, A=9), this effective nucleon being defined by the combination '4/9 
proton -t- 5/9 neutron'. The theoretical calculations use the CTEQ5M set of parton distribution 
functions (pdf) for the proton pdf's, and assume that the neutron pdf's are simply deduced from 
the ones of the proton by isospin symmetry. The KKP set [Q is used for the fragmentation functions 
into a pion. Our theoretical results are given for a common scale choice for the renormalisation 
scale ^, initial state factorisation scale M and final state fragmentation scale Mf equal to n = M = 
Mf = (3/4)/!" where /i° — {Pti + Pt2)/2 is the average value of the pt's of the two pions. At 
fixed target energies, the scale uncertainties remain pretty large in NLO calculations. We made a 
rapid investigation]^ of the scale dependence of the predictions around the above choice, using also 
the two other followings sets of scales: /i = M — (l/2)/i°, Mf — and ^ — M — Mf — jjP. 
The choice ji = M = Mf = (3/4)/i° turns out to accommodate to the data better in the case of 

^The inclusive hadroproduction of pion pairs has also been studied earlier by the CCOR ISR experiment |l7| and 
the NA24 fixed target experiment at CERN. We choose to focus on the data collected by the E706 experiment, 
which is more recent and gathers higher statistics. 

comprehensive procedure would be, as in one particle inclusive production, to perform a detailed study of these 
dependences with respect to /i, M , Mf separately. This would help, at least, in quantifying the scale uncertainties 
more quantitatively. This would also possibly detemiine a region of local stability with select a set of optimal scales, 
invoking the Principle of Minimal sensitivity [pcj, |2l|. At fixed target energies, such an optimal scale choice can be 
found in the case of sinp- ln p ion inclusive case. However, the surrounding region of local stability is a narrow area 
sitting on a sharp ridge |22[ , an unsatisfactory situation. Furthermore, this type of study is quite computer time 
consuming already in the single particle inclusive case. Its cost becomes prohibitive in the case of correlations in 
dipion production. Therefore we did not perform such a complete study. 
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the pt distribution of each pion, as can be seen in Fig. whereas the higher (loweiQ) scale choice 
undershoots (overshoots) the data, cf. Fig. |l|. In addition, we made comparisons with E706 data in 
pp collisions with an energy beam of 800 GeV (v^ ~ 38.76 GeV), using the same scale choice and 
reaching similar conclusions. This somewhat legitimates the choice of scales adopted. In this section 
we display the comparisons between the data and our NLO estimates for the various observables 
measured by the E706 experiment. A corresponding phcnomcnological analysis will be presented in 
the next section. We see that these observables split qualitatively into two groups, according to the 
degree of agreement or disagreement with the experimental spectra. 

3.1 Well described observables 

For the given scale choice, the agreement is reasonably good in shape, as well as in magnitude, for 
the distribution of transverse momentum of each pion in pairsj] da/dpT (Fig. |^), for the invariant 
mass distribution of pairs da /drriT^T^ (Fig- ^ - except in the very first bins at low invariant mass -, 
for the rapidity distribution of the pairs da/dyT^^ (Fig- §), and for the distribution in cos6'* (Fig. 
The variable cos 6* is defined by^ cos6'* ~ tanhy* where y* = (j/i — j/2)/2 is half the difference of 
the rapidities of the two pions in a pair. 

The reasonable agreement observed was expected for these observables over the whole measured 
spectrum. Let us however focus in particular on the distribution da/dpr- The remarkable fact is that 
the agreement is fairly good over the totality of the measured spectrum. This may be contrasted with 
the situation in the case of one-particle inclusive production, measured also by the same experiment, 
and which was confronted to a NLO calculation in Ref. Let us first briefly remind the results 
which were obtained in Ref. . The calculation used the BKK set of parton to pion fragmentation 
functions ||^. The scales, selected by the application of the Principle of Minimal Sensitivity ||20|,pl[|, 
were close to Pt/3- This set of scale somehow maximized the theoretical prediction, and turned 
out to minimize the discrepancy between the latter and the data. Above pr > 5.5 GeV, no large 
difference between the shapes of the NLO calculation and the data was found. On the other hand, 
the NLO calculation undershot the data by about 50%. The discrepancy between the NLO result 
and the data was found larger with other, especially larger scales (a factor 1.8 for the scale choice 
Pt/2). Furthermore, in the lower range of the pr spectrum below 5.5 GeV, the slope of the data is 
increasingly sharper than the one of the NLO post-diction as the pr decreases. On the contrary in 
the case of the pr distribution of each pion in pair production, the agreement is reasonably good in 
shape, and also in magnitude with the intermediate scale choice fj. = M = Mf = (3/8) (Pti + Pt2) 
used. One notices that the above scale choice involves higher values (by nearly a factor 2) than the 
optimal scales in the case of single inclusive production. The agreement in shape is rather robust 

^One cannot go down to too low scales. Another choice of scale commonly used would he /i = M = Mf = 
(PTi -l-PT2)/4 (half the average value of the pt's of the two pions). However the lowest factorisation scale Mq at 
which the KKP set starts is Mq = \/2 GeV. With the selection cut pTmin = 2.5 GeV on each pion, this choice 
would lead to a scale that could fall below Mq. Technically, a backward evolution, or, more drastically, the freezing 
of the pdf 's below Mq can be enforced. However the use of such questionable remedies would induce extra sources of 
theoretical uncertainties. 

'^According to this definition each event enters twice, i.e.: 

d(T d(T da 
dpT dpTl dpT2 

®This variable has been initially introduced in the literature in the case of 2 ^ 2 processes, where d* is the angle, 
defined in the center of mass system of the pair, between the direction of the beam and the direction of the pions, 
emitted back-to- back in this frame. This definition is convenient at lowest order (LO). It is however difficult to extend 
it in a simple and unambiguous way to processes with more than two particles in the final state, as it is the case 
in calculation beyond LO in particular. In the case of a 2 — > 2 process, cos 9* turns out to be equal to tanhy*. 
In addition to provide an alternative definition which is manifestly invariant under longitudinal boosts, this identity 
allows to extend the definition of cos 9* to processes with more than two particles in the final state. This is the 
definition used by the E706 experiment [hi, and the one which we adopt here. 
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with respect to variations in the arbitrary scales /i, M and Mj, at least within the three common 
values which we scanned, cf. Fig. |l|. Furthermore, the agreement in the case of pair production 
holds even down to fairly low values of pT close to pTmin — 2.5 GeV above which the events 
were experimentally selected, in sharp contrast with the situation in the single inclusive case. This 
fact may appear somewhat surprising, and even counterintuitive because one may think that the 
agreement should be better, the more inclusive the observable. This puzzle will be examined in some 
details in the next section. 

3.2 Poorly described observables 

In contrary, the NLO approximation is completely inadequate to account for the distribution of az- 
imuthal angle difference between the two pions, cJa/d^^^ (Fig- H), and for the distribution da/dpT out 
(Fig. 1^. The variable PTout is defined as follows. The direction of the beam axis and the transverse 
momentum of one of the pions define a plane; the component of the px of the other pion which 
is orthogonal to this plane defines the PTout of this pioriPl The disagreement is also sharp for the 
transverse momentum distribution of pairs da/dq^ (Fig. The transverse momentum balance dis- 
tribution da/dZ (Fig. ^ is not adequately described by the NLO approximation either, particularly 
in the region Z < 1. The px balance variable Z is defined for the pion j by|^ Z — — (pVi • PT2 ) / \PTj \ ^ ■ 
The shapes of these observables, especially of ^tttt and pTout distributions are notoriously infrared 
sensitive respectively when (pT^Ti tt and pTout 0. Therefore it is not surprising that the NLO 
calculations do not describe these regions adequately. The noticeable fact is that even the tails of 
these distributions are unmatched. The NLO estimates fall much too rapidly compared with the 
data, the former being more than one order of magnitude below the latter. These features will be 
examined in detail in the next section. Last, in Fig. we notice also the poor account of the very 
first bins at low invariant mass in the invariant mass distribution by the NLO calculation. 

4 A classification of dipion observables 

The observables of inclusive pion pair production can be theoretically classified into two main cate- 
gories, according to their respective infrared sensitivity. This theoretical classification matches with 
the phenomenological one above, according to the 'degree of agreement (or disagreement)' between 
the NLO predictions provided by DIPHOX for inclusive neutral pion pair production and the recent 
high statistics data of the fixed target E706 experiment at Fermilab. 

The first category which we will consider gathers the observables which are infrared insensitive. 
By this we mean that the contribution associated with the multiple soft gluon emission is not 
logarithmically enhanced order by order in perturbation theory. Typical representatives of this class 
of observables are: the distribution da/dpx', and da/dniTrTr, at least when the invariant mass m,rir 
is larger than 2pTmin (the reason for this restriction will be examined later). Other examples are: 
the distribution da/dcosO*; and the rapidity distribution of the pairs, da/dY. 

The opposite class gathers the infrared sensitive distributions. The shapes of these distributions are 
well known to be strongly influenced by multiple soft gluon emission in the vicinity of some 'critical 
point'. Depending on the observable considered, this critical point may lie at the elastic boundary of 
the spectrum; or it can stand inside the physical domain. Order by order in perturbation theory the 
infrared sensitivity is signed by the appearance of large logarithms of Sudakov type which become 
infinite at the critical point. Among these distributions, one may distinguish between two types. 
A first family is made of the distributions for which the perturbative calculation is singular at the 

^In other words, the Pt out of pion j = 1,2 is equal to p^jsint/i^^. Notice that each event enters twice in this 
distribution. 

^"^In this observable as well each event enters twice. 
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critical point order by order in perturbation theory. In the present case of pion pair production, 
typical representatives are da/d(j)T„y, da / dpT out- Any fixed order calculation is certainly inadequate 
to estimate this type of observables in the vicinity of the critical point, and, whenever possible, a 
Sudakov type resummation has to be carried out to restore a predictive power of the perturbative 



QCD calculation. As will be commented in subsect. 4.2, the colour structure of the hard subprocess 
in hadroproduction of pairs of hadrons is more complicated than in the celebrated case of the Drell- 
Yan process. This makes a soft gluon summation beyond the leading logarithmic accuracy quite 
sophisticated in this case. We do not provide a study of these observables based on a resummed 
calculation. A second family gathers observables which are less infrared sensitive than the previous 
ones, in the sense that their fixed order estimates do not diverge at the critical point; therefore the 
inadequacy of fixed order calculation to account for such observables is more difficult to foresee. 
Among these observables, let us mention the distribution of transverse momentum of pairs, da/dqx, 



and the distribution da/dZ where the variable Z has been defined in subsect. 3.2 . They would belong 
to the first family, if the fragmenting partons in the final state were, instead, observed particles - 
as, e.g. it would be the case if they were prompt photons. In the collinear fragmentation model 
used here, the screening of the infrareed sensitivity stems from the convolution with fragmentation 
functions, this provides extra integrations over longitudinal fragmentation variables. 

Computing the NLO truncation for infrared sensitive observables is more than an academic exer- 
cise. Indeed, whereas the vicinity of the critical point is infrared sensitive, the tail of the distribution 
may be expected safer for fixed order perturbative calculations. Yet the infrared sensitivity of these 
observables connects the perturbative and non perturbative regimes. At fixed target energies, de- 
pending on the observable, on the selection cuts, etc., the non perturbative effects may happen to be 
not confined to the vicinity of the critical point, but instead can spread over a large fraction of the 
physical spectrum, if not the whole. A confrontation of the NLO truncation to the measured spec- 
trum can give a hint of the depth to which infrared, in particular non perturbative, effects penetrate 
a distribution away from the infrared critical point. This is what we have in mind while performing 
this confrontation here. Another interest concerns observables which are infrared insensitive apart 
from a few bins where, for example, kinematic cuts make them locally semi inclusive. The discrep- 
ancy between the NLO prediction and the measured spectrum in these bins can be related to some 
infrared sensitive distribution inadequately accounted for by fixed order truncation. This is typically 
the case with the lowest bins of the invariant mass distribution, where the discrepancy between NLO 
and data can be understood in the light of da /dcj}-,^-^, as will be explained below. 

Let us add a comment about the classification which we present here, especially regarding the infrared 
sensitive observables. Phenomenological models are sometimes used for 'effective' kx dependent 
distribution - and fragmentation - functions. For example, a popular procedure 23| revisited 



recently [ p4[ assumes a x independent, gaussian kx distribution functions for on shell partons. In 
such a model, the behaviour of distributions such as da /dipTr-K and da/dpTout at respectively (/),r7r tt 
and pt out are smeared by the convolutions with the kr distributions. Therefore, according to 
such a model, all distributions would be affected by kx effect, i. e. the first family would be empty. 
In order to avoid any confusion in the terminology used, let us emphasize the following point. The 
models of kx smearing mentioned above depart from the property of collinear factorisation, and they 
do not account for the higher orders corrections to the hard subprocess in perturbation theory in a 
consistent way. On the contrary, the classification which we present here follows closely the property 
of collinear factorisation which holds in perturbative QCD. 

4.1 Infrared insensitive observables 

Transverse momentum: one-particle inclusive vs. inclusive pair production 

Let us analyse the difference between the px distribution of each pion in the present case of dipion 
production, and the px distribution in one particle inclusive production studied in Ref. [E2|. We 
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have to disentangle two issues. The one concerns the difference in the ingredients used in Rcf. 
vs. the present work. It will affect the whole spectrum, and it is especially relevant in the larger pT 
range. The second concerns the physics involved, and will particularly affect the lower range of the 
spectrum. We consider both of them successively. 

First, the sets of fragmentation functions used in either case are different. This difference turns out 
to be quite significant in the case of single pion inclusive production. It affects the whole spectrum, 
and it matters especially in the higher range of the px spectrum, pr > 5 GeV, where we expect 
perturbative QCD to apply. Indeed, as seen in Fig. the u and d fragmentation functions of the 
KKP set are larger than the corresponding ones of the BKK set: more than 50 % above z = 0.7, and 
the ratios increase as z increases towards 1. The difference (in the same direction) is even larger for 
the gluon fragmentation: more than a factor 2. This is possible because the sets of fragmentation 
functions are extracted from data coming from e^e^ annihilation into hadrons which constrain 
the fragmentation functions essentially in the range 0.1 < z < 0.7 and because, unlike pdf's, the 
fragmentation functions are not tightly constrained by a momentum sum ruleQ. 

It is worth noticing that for pr ranging from 4 to 10 GeV, the average value < z > of the fragmen- 
tation variable z in the one particle inclusive case is < z >~ 0.75 — 0.90. Therefore the theoretical 
estimates for single inclusive production rely on extrapolations of the fragmentation functions out- 
side the region where they are actually constrained by the data from which they are extracted. 
Consequently, as emphasized in Ref. [g^, the NLO predictions are plagued by a rather large un- 
certainty. The effect of the replacement of the BKK set by the KKP one is shown in Fig. ^ The 
discrepancy in normalisation between the theoretical prediction and the data is strongly reduced, 
for the scale choice pt/3. 

Comparatively, the theoretical calculation of inclusive pion pair production is less affected by the un- 
certainties on fragmentation functions at large z. The reason is that the values of the fragmentation 
variables zi , Z2 which dominate the production of pion pairs extends in a region of smaller values 
than the < z > in the single inclusive case, zi^2 typically ~ 0.5, where the e'^e~ data constrain the 
fits of the fragmentation functions quite well. This feature is illustrated by Fig. This can be 
traced back to the fact that the region in z, which dominates, results from a competition between the 
cost in energy momentum involved in the partonic subprocess - the pdf's tend to favour the smallest 
possible X fractions for incoming energy-momenta and therefore larger values of z - and the price 
to pay for a fragmentation near the end point z = 1 - the fragmentation stage favours small values 
of z. The production of a pair of two hard particles involves two fragmentation functions instead 
of one as in the single inclusive case, which enforces the role of fragmentation in the competition. 
The situation will be very different regarding the reducible background to Higgs boson search at 



LHC 1 38 , because isolation requirements will select the large z region. This issue will be addressed 
in detail in a forthcoming article [p^ . 

The difference in the dominant ranges of z in one particle inclusive vs. pair production sheds also 
some light on the apparent difference in scales which seem to be preferred by the data in either 
case, the optimal scale choice ~ Pt/3 in the one particle is apparently smaller than the choice 
(3/8) {pti + PT2) used in pair production. Let us consider the inclusive production of a single 
pion in the region where the E706 data agree reasonably well with the NLO calculation using the 
KKP set, say, near px = & GeV. Such pions probe a partonic regime with a transverse momentum 

^'^In the framework of independent fragmentation functions, there exists indeed a sum rule; 



J2 I' dzzDH/a{z,Mj) : 



The latter means that a parton of given species a fragments into any hadron species H with probability 1. It is 
however much less constraining in practice than the one holding for pdf's which relates the pdf's of all the species for 
a given incoming hadron. 
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~Pt/<z>^7 GeV, and the optimal scale is ^ 2 GeV. A pair of pions, dominantly back to back, 
with individual ^ 3 GeV probes a parton regime with parton transverse momenta ~ pt/ < ^1.2 > 
comparable with the previous one; we notice also that the scale choice used in this case is 9/4 ~ 2 
GeV. Consequently, the same partonic kinematic regime is probed in these two cases, and are found to 
be consistently in reasonable agreement with the same theoretical framework. We therefore conclude 
that the situation is actually comparable between the two cases in this pr range. Of course this 
rapid investigation considers only the comparison with the data of E706. To obtain a comprehensive 
understanding, one needs to perform a complete study revisiting single inclusive production at fixed 
target energy involving other fixed targets experiments. 

Let us now examine the difference in the behaviours of the pT spectra of single inclusive production 
vs. each pion in pair production in the lower range of the pr spectrum. The NLO calculations 
performed in the QCD improved parton model according to the property of coUinear factorisation 
neglect non perturbative, power corrections ~ O (1 GcV/pt)") where n depends on the observable 
considered. Therefore this 'leading twist' approximation is valid when pT ^ 0(1 GeV). On the 
other hand, NLO calculations are not expected to predict correctly inclusive pt spectra when p^ 
becomes too small, and distortions of the perturbative QCD predictions are expected up to ~ a 
few GeV. One may instead wonder why the NLO calculation and the data agree well on the slope 
of the distribution of each pion's pt for pair production down to pT ''^ PTmin = 2.5 GeV. 

Among the neglected effects, the ones induced by finite transverse momenta (kx) of incoming partons 
are often assumed to be the most important. This concerns in particular the understanding of the 
Pt distribution at low pr in one particle inclusive production due to the so-called 'trigger bias' effect, 
tentatively visualized in Fig. A current practice models the latter by assuming an effective kr 
distribution function]^ for the incoming partons, still considered on shell, and convolutes this with a 
LO calculation. If one notes g{kx) this kx distribution, and assuming that it depends only on IfcT^I, 
one has schematically 

^(fkTgikr) 



dpT J \pT - h 



\ 2 pi ^""y p% 

This means that the contribution from subprocesses initiated by partons with incoming transverse 
momentum pointing in the direction of the measured pion is enhanced, and this leads to an increase 
of the cross section in average, especially at low pT. Let us however stress that this model is not 
consistent^ with the property of coUinear factorisation, established in perturbative QCD in the 
leading twist approximation, for inclusive processes with hard probes. Attempts of a consistent 
formulation have been proposed recently |2^,^^. This formalism is still in an exploratory stage, and 
no firm and quantitative conclusions can be drawn yet. 

On the other hand, the distribution of each pion's pt in pair production involves the detection of two 
hard particles, dominantly back-to-back in the transverse plane. Therefore, any fc^ kick which would 

^^For example, the popular gaussian model mentioned above. 

^^In addition it requires the use of several ad hoc extra parameters and cuts, in particular to prevent the convolution 
integral in eqn. (hi) from diverging. The divergence of eqn. (hi) when fcy = Pt corresponds to a soft partonic subprocess 
boosted transversally by the fey kick. The divergence obtained is an artefact of the model. Another cut-off has to be 
introduced to prevent the longitudinal momentum of an incoming parton a and the momentum of its parent hadron 
H from being in opposite directions. This nasty can occur at non zero fcy since 

^ i 4„ 

where Xa is the light cone momentum fraction of H carried by a. This problem is an artefact of the infljdelling (an on 
shell parton with non zero transverse momentum and a Lorentz invariant momentum fraction). See |23[ for a review. 
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favour one pion in a pair would penalize the other Q . This is tentatively pictured by the top drawing 
in Fig. In particular, in the lower range of the px spectrum, the pT of the penalized pion may 
even be pulled down below the selection cut pTmin so that the event may be discarded, cf. bottom 
drawing in Fig. Consequently, when the selection cut is symmetric, the mentionned trigger bias 
is reduced. If the mechanism of kx gives the dominant power correction in the single inclusive px 
spectrum at low px, its lesser influence on the px spectrum of each pion in pair production allows 
perturbative QCD to extend its range of validity down to fairly low px values. Let us notice that 
the absence of trigger bias effects is expected as well if, instead of the distribution of each pion's 
Px, but with the same selection cuts, one would consider the px distribution of the 'leading' pion 
(the one with the largest px in the pair), even though the latter seemingly displays the px of only 
one particle. On the other hand, an interpolating case would be the distribution of each pion's px 
in pair production when the selection cuts on the two pions are asymmetric. The behaviour of the 
distribution of each pion's px is expected to exhibit the same trigger bias effect as the single inclusive 
case when the px cut on the second pion becomes very low. 



The invariant mass distribution of pairs, da/draT^T^ 

The D/PiJOX calculation and the E706 data agree reasonably well on this observable for mTrn- > 5.0 
GeV. On the other hand we see a sharp disagreement on the first two bins in invariant mass. In this 
respect we notice that 



m^TT = V2pTiPT2(cosh(2y*) - cos0^^) (1) 

(pT^TT being the azimuthal angle between the two pions. At LO, the 2^2 kinematics forces the px 
of the two pions to be back-to-back (cos^^jr = — 1) and the value 5 GeV, which corresponds to 2 
PTmin, IS the minimal value which can be taken by mTr-n- at LO. In other words, the first two bins 
of the theoretical calculation are filled only by NLO, inelastic contributions, for which </i7r7r < tt. 
The sharp disagreement between the NLO calculation and the data in these two bins is therefore 
correlated with the situation regarding the distribution in azimuthal angle difference between the 



two pions in a pair, which will be examined in sect. 4.2 



The distribution da /d cos 9* 

As seen on Fig. |^, the measured distribution da /d cos 0* and the corresponding NLO distribution 
show a fairly good agreement in shape but an overall discrepancy in normalisation|^, the theoretical 
prediction being systematically above the data by a constant multiplicative factor ~ 1.4 for the scale 
choice used. One may wonder why, for the given scale choice in the NLO calculation, and for the 
same data sample, the agreement in normalisation seems somewhat worse for this distribution than 
what is found in average for the other infrared insensitive distributions. 

Firstly, it should be stressed that this cos 6* distribution is much more tricky than its naming 
suggests. Indeed, the distribution shown in Fig. || is not literally the bare quantity da /d cos 9*. Let 
us recall that the 'naive, bare' distribution da/dcos9* , without additional cut, has a bell-like shape 
peaked at y* = and decreasing towards as |j/*| is increasing. This results from a competition 
between the behaviour of the partonic matrix element squared, and the fall-off of the pdf's when 
X ^ 1. The matrix element squared behaves as as (1 ± cos6'*)~" when cosO* increases, in the case 
of a massless spin s exchange in the i or u channel: a distribution unveiling this behaviour would 
namely probe the short distance dynamics. On the other hand Xi ^ 1 for i = 1 or 2 when \y*\ 
increase^, and this effect wins. Additional cuts are implemented to freeze the influence of the pdf's 

A current presentation of this observable normalises the measured distribution and the NLO prediction separately 
in such a way that the bin content corresponding to cos 6* = is equal to 1. The discrepancy in normalisation is 
therefore hidden by this procedure. 

^^See for instance psj for explicit LO expressions. 
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and magnify the short distance dynamics. The additional cuts restrict the data sample providing 
this cos 9* distribution: this sample is not the same as the one yielding the other infrared insensitive 
distributions, but only a subsample of it. In particular, an auxiliary selection cut on the invariant 
mass m^TT > 7 GeV is imposed, in order to avoid pT-threshold effects induced by the selection 
cut pt > PTmin ■ The whole distribution carved out this way is dominated by the pairs with 
m^TT ~ 7.0 — 8.0 GeV. We then notice on Fig. for the invariant mass distribution that the ratio 
data/NLO in the corresponding bins is approximately 0.7, a value consistent with what is observed 
on Fig. ^ bottom for the cos0* distribution. In summary the normalisation of the whole cos0* 
distribution is given by the one of only a few bins of the m^jr distribution. This makes the overall 
normalisation in the NLO cos 9* distribution very much scale dependent. If, instead of choosing 
the scale 'preferred' by the each pion's px distribution in the theoretical calculation as explained in 
sect. H, we take a scale which makes the ratio Data/Theory close to 1 for the m-TrTr distribution, this 
apparent normalisation problem disappears. 

4.2 Infrared sensitive observables 

We now successively examine the two families of infrared sensitive observables which we identified at 
the beginning of this section. We first focus on the distribution dcr/cZ^^Tr as an example of observables 
with strong infrared sensitivity. Next we consider the cases of the qt and Z distributions, as examples 
of observables for which the infrared sensitivity is smeared. Later we examine how the sensitivity of 
these types of observables to non perturbative effects evolves with increasing individual transverse 
momenta of the pions and/or with increasing invariant masses of the pairs. 

The angular distribution dcr/d^Trir 

As can be seen on Fig. the 'NLO distribution' and the experimental spectrum are absolutely 
different. The LO contribution to the 'NLO distribution' is proportional to S{Tr — (pT^Tr), and the 
inelastic part of its NLO contribution diverges as ©(ofg lii(7r — (j)^^)/[T: — (/itttt)) when (/itttt — > tt. 
Therefore the 'NLO distribution' is mainly concentrated in the last bin containing (pT^-j^ = tt at the 
elastic boundary of the spectrum. In contrast the experimental distribution is spread out much more 
and flatter. The discrepancy is so spectacular that one has to ask the question, why there is such 
a satisfactory agreement for the pT distribution, cf. Fig. H, and, similarly what the reason is for 
the sharp disagreement in the first two bins for the invariant mass distribution in Fig. |^, while the 
agreement is satisfactory in the other bins. 

We notice that, although the shapes of the 'NLO distribution' and of the measured distribution 
are completely different, the area below these distributions are approximately equal. This suggests 
that long distance phenomena, which are ignored in the NLO calculation, proceed to nothing but a 
drastic kinematic 'reshuffle' of the configurations of the hard partons involved in the subprocess into 
the ones of the measured pions, i.e. a probabilistic redistribution of kinematic variables occurring 
with total probability 1. This hypothesis is supported by the following argument. Looking back to 
the invariant mass distribution of pairs, we noticed previously a strong disagreement between the 
NLO prediction and the data in the first two bins. Quantitatively, the theoretical result is more than 
a factor 25 smaller than the measured central value for the second bin 4.5 GeV < ra-^T^ < 5.0 GeV, 
and it is about a factor 60 smaller in the first bin 4.0 GeV < TOtttt < 4.5 GeV. Given the fast fall-off 
of the Pt distribution the pion pairs for which the individual pions' pt are minimal i.e. close to 2.5 
GeV both dominate the first bins of the invariant mass distribution, and control the (/)^7r distribution. 
Furthermore, pairs with small y* dominate. Therefore, approximating the identity by roughly 
m^^ 2prmi„ sin((/)^7r/2) we infer that the first and second bins in invariant mass are dominated 
by pairs with respectively ^tttt ~ 115 ° and (pTr-w ~ 140 °. If we now look at the discrepancy between 
the 'NLO calculation' and the measured distribution for the ^tttt spectrum at these values of (j^TrTv, 
we find respectively ~ 40 and 20. These values are comparable with the ones for the discrepancies 
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between NLO results and data for the invariant mass spectrum in the corresponding first and second 
bin in invariant mass. 



A proper understanding of the shape of the distribution da/dcj)^^^^ requires the all order summation of 
the large Sudakov logarithms which appear order by order in the theoretical calculation. The latter 
come from soft gluon emission in both initial and final states - and, beyond leading logarithmic 
accuracy, from interferences between them. The colour structure involved |2^,|2^ is quite more 
intricate than in the well-known processes of Drell-Yan [2g|-|3^ and nearly back-to-back hadron 
production in hadronic e'^e~ annihilation ||3^ . This summation has not been performed yet for the 
(j)TnT distribution. Beyond this summation, non perturbative effects are also involved. In the analysis 
of infrared sensitive distributions of electrodynamic probes in fixed targets experiments, such as 
the transverse momentum distribution of lepton pairs in Drell-Yan or of photon pairs in prompt 
photon pairs production js^ , such non perturbative effects have been shown to be important over 
a large fraction of the spectrum. In the present case they may have even larger effects. They may 
come in particular from some kx kick on initial partons, as well as on the hadrons produced by the 
fragmentation process. Since both the initial and final states are affected by these phenomena, and 
since the selection cut on the px of the pions is not much larger than the typical scale ~ 0(1 GeV) of 
non perturbative phenomena in strong interaction, even a rather moderate effective fc^ kick - several 
hundreds of MeV 'per hard coloured leg' - may be able to reshuffle the hard parton configuration 
drastically, as one can convince oneself with the cartoon drawn in Fig. |l^. 

It is important to understand how it is conceivable that the infrared insensitive observables can be 
reasonably estimated using perturbative, fixed order calculations, whereas there exist distributions, 
such as da/dipT^T^ which are so heavily distorted by non perturbative effects. The key point is the 
following: if this probabilistic reshuffle acts on kinematic variables which are integrated over during 
the calculation of some distribution, and if the available phase space for this kinematic variable is not 
severely amputated by some selection cut in the range considered, this reshuffling leaves nearly no 
fingerprint in this observable. This seems to be what happens for the distribution of each photon's 
Pt (the variable cjiTrn is integrated out when the pT of the other photon is integrated out). The same 
conclusion holds also, for example, for the invariant mass distribution as soon as the invariant mass 
is large enough compared to 2pTmin- It does not hold when rriTrTr is near or below this value because 
of the restriction acting on (fi-^^^ through eqn.(^. 

Transverse momentum distribution of pairs, da/dqr 

As we mentioned in the general introduction of this section, the situation in the case at hand is 
different from the one usually faced concerning this observable in the Drell-Yan process/production 
of a weak vector boson, or in the production of pairs of prompt photons]^. E.g., in the Drell-Yan 
case, the calculation of this distribution is singular at = order by order in perturbation theory. 
In particular the zeroth order contribution to dcr/dq^ is proportional to S{q^) by transverse momen- 
tum conservation, and logarithmically enhanced '-I-' distributions are involved at every successive 
order. On the contrary, in the present case, the double convolution with the two fragmentation 
functions makes the calculation finite at every transverse momentum qt of the pair, order by order 
in perturbation theory. In particular it is non vanishing already at lowest order for every qt > 
allowed by kinematics. One may therefore expect a somewhat reduced infrared sensitivity, compared 
with the (/)^7r distribution. The comparison between the NLO calculation and the measured qt dis- 
tribution, cf. Fig. ^, namely the examination of the discrepancy between the former and the latter 
suggest that it is hardly the case. Whereas the discrepancy reaches quantitatively a magnitude of 
20 to 50 in the ^tttt case, it is 'only' about a factor 5 to 10 for the qx distribution. Moreover the 

More precisely, we mean here the production of photon pairs through the so-called 'direct' mechanism, cf. sect. |2| 
and Ref. Iisl. On the contrary the production of pairs of 'photons from fragmentation' is very similar to the present 
case. 
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NLO calculation is still inadequate in shape as well over all the measured qr- This inadequacy is 
related to the discrepancy in the (/)^7r distribution examined in the previous subsection. Let us write 
~ {pti ~ PT2)^ + 2pTiPT2(l + cos0^^). The cjiT^^^ distribution is concentrated in the vicinity of 
(j)^^ = TT whereas the experimental one is more flat and spread. Therefore the NLO result for the qt 
is peaked near qt = and is higher than the experimental spectrum in this region. As for configu- 
rations with a larger qt, they can have two origins. Either the individual pt's have to be larger, but 
NLO and data agree on the px distribution, hence this origin contributes by the same amount to the 
NLO result and to the data; or (/)^7r is far from tt, in which case the NLO 0^7r distribution is much 
below the experimental one. This makes the NLO qt distribution fall faster than the measured one 
as qr gets larger. 

Transverse momentum balance distribution da/dZ 

This observable is another case of distribution which is singular e.g. in diphoton production, namely 
at the critical point Z = 1. Indeed, in diphoton production, the contribution of the 'direct' mecha- 
nism is proportional to i5(l — Z) at lowest order by transverse momentum conservation, and involves 
logarithmic singularities at every higher order. Due to the same reason as in the previous exam- 
ple, this distribution is, instead, finite for all Z, order by order in perturbation theory. As can 
be seen on Fig. ^, the comparison between the NLO calculation and the experimental distribution 
roughly agrees for Z > 1, but is inadequate for Z < 1. This can be understood as follows. From 
Z = — pti/pt2 cos 0^7r one sees that pxi > pt2 is a necessary condition to have Z > I. Since this 
occurs preferably with values of 0,r7r not far from tt, one understands that the NLO curve overshoots 
the data slightly regarding the last bins of the distribution, cf. Fig. |^. Yet the reversal of the 
ordering of the respective magnitudes of the NLO and measured (jj-^Tr distributions occurs for a value 
4>TrTT close to TT. The integration over the fragmentation variables averages these two opposite or- 
derings, and makes the discrepancy between the NLO and measured Z distribution rather small for 
Z > 1. On the other hand, for Z < I, pxi and pT2 are unconstrained and most likely of comparable 
size. This Z region thus selects the lower range of cos ^tttt and for small Z one finds a similar large 
discrepancy as for the (^,r7r distribution in the lower (j)Tnr range. 

4.3 Energy dependence of the non perturbative effects 

The sensitivity of the infrared sensitive observables examined above to non perturbative aspects is 
strong and extended to a large fraction, if not the totality, of the measured spectrum. This is mainly 
because the E706 data were selected above a cut pTmin > 2.5 GeV, and because the normalisations 
and shapes of these distributions are dominated by the pairs in which the individual pt's of the pions 
are close to pTmim a value not large with respect to the typical scale of non perturbative effects. 
As we are interested in making predictions for collider experiments, especially for the LHC ||l2|, it 
would be interesting to see how the situation evolves with the scales. Two situations have to be 
distinguished: 1) the increase of the 'hard scale', loosely speaking - i.e. the invariant mass of the 
pair, or the individual pr's of the pions in pairs - at fixed V^, and 2) the increase of at fixed 
'hard scale'. 

In the archetype of the Drell-Yan process, for the distribution da/dqTdm of lepton pairs, this question 
has been studied relying on resummed calculations. In particular, the authors of Refs. |^,^ have 
studied how the sensitivity to non perturbative effects at gT = evolves as a function of m. They 
determined which region dominates the Sudakov-type factor in impact parameter (b) space, by 
mean of a saddle point analysis. They showed that the dominant region in b space is purely non 
perturbative at m values in the fixed target energy range, but that it migrates towards smaller values 
of b and ultimately fits completely inside the perturbative region Kqcob ^ 1 when m becomes large 
enough. This indicates that the non perturbative effects become confined to a closer and closer 
vicinity of the critical point as the invariant mass m increases. This feature makes it difficult to 
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extract a precise parametrisation of these non perturbative effects from experimental data 
On the other hand it imphes that theoretical predictions at higher, coUider energies are less affected. 

In this subsection, we perform a phenomenological study through which we try to extract informa- 
tion about item 1) within the range covered by the E706 experiment. As noticed in the last two 
subsections, the ^tttt distribution is the key of the understanding of the others. A useful information 
could be extracted from data samples ordered according to increasing values of px mm for fixed v^. 
Tabulations of data of this type are available for the invariant mass distribution. However not much 
information can be learned from it since this observable is blind to non perturbative effects over 
most of its spectrum. On the other hand, such tabulations are not available for the more instructive 
07r7r distribution. An alternative could be provided by a refined double binning of the data according 
to the infrared sensitive variable considered and, simultaneously, to increasing invariant mass, the 
latter corresponding loosely to increasing individual pt's of each pion in pairs. This type of double 
binning is available for the pTout distribution, but again not for the (/)^7r distribution. The only 
information available about (/)7r7r in this respect is the average < (/)^7r > • 

Increasing the hard scale at fixed \/S. 

Let us thus consider the average values of the azimuthal angle difference between the two pions 
of pairs, < ^^^r > as a function of the invariant mass of the pair. Of course the shape of the (/jt^t^ 
distribution is infrared sensitive near 0^7r = and, to any fixed order in perturbation theory, the 
theoretical estimate of the ^^^r distribution is divergent at (/)^7r — i^- However when the average 
< (p^^ > is computed, the Sudakov logarithms are smeared out by the integration order by order in 
perturbation theory. There is no logarithmic enhancement in < (p-^T^ > order by order in perturbation 
theory anymore. Consequently, a fixed order calculation of < (p-m, > is expected to become more and 
more accurate for harder and harder pion pairs, without a detailed knowledge of the shape. Non 
perturbative effects are expected to induce deviations in the form of power corrections O (1 GeV/pr), 
©(IGeV/m^^). 

For m^r^ > 5 GeV, the full range 105° < 4> < 180° is accessible]^, and < (p^^T^ > is dominated 
by the nearly back-to-back region. This dominance of the region ^tttt ~ tt is actually stronger in 
the NLO result than in the data. Indeed, the pairs of pions with low individual px and relatively 



low invariant mass dominate the </)7r7r distribution discussed in subsect. L2, The latter is thus 
strongly sensitive to non perturbative effects, which make it much flatter than the shape of the NLO 
result. Therefore the experimental values of < > should be lower than the NLO values. The 
distribution da / dm-^T^dcj)-,^-^ should then become steeper and more concentrated near the back-to-back 
region with increasing invariant mass, if non perturbative effects become less and less important. 
The experimental values of < (p^^^ > should increase with increasing invariant mass, and the difference 
between them and the corresponding NLO results should decrease. The general trend of the data 
seem to agree with this expectation, as suggested by Fig. |l^. To be fair however, the very large 
error bars prevent from drawing any firm conclusion. More generally, the first moment of infrared 
sensitive distribution is not sensitive enough to long distance effects. Cumulants of higher orders 
of these distributions would be much more instructive in this respect, not speaking about these 
distributions themselves. 

We examined also the distribution da / dpr outd'mTT-K sliced in broad bins of invariant mass. As seen 
on Fig. ^ all plots for each separate bin of invariant mass are similar to the one shown in Fig. 
The NLO result is inadequate to account for the measured distribution for all invariant masses, and 
nothing instructive can be extracted from it regarding the range of influence of non perturbative 

^^The situation in the first bins with m-n-K < 5 GeV is pecuhar. In these bins the selection cut pT > PTmin forbids 
the back-to-back region as discussed previously. The average < (j>^^ > is thus much smaller than n. This holds both 
for the NLO result and the data, and no relevant information can be extracted regarding the problem concerned by 
this discussion. 
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effects with respect to the hwariant mass. Last, in the case of the distribution da/dZdrnj^j^, we 
notice on Fig. ^ that the discrepancy between the NLO result and the data for Z < 1 decreases 
with increasing TOtttt- This holds at least for Z values not very small with respect to 1, for bins in 
invariant mass which gather enough statistics so that any conclusion can be drawn. 

It would have been interesting also to identify how the magnitude of non perturbative effects evolves 
when the beam energy increases from 530 to 800 GeV. Unfortunately, the level arm in a/S is too 
small and the errors bars are too large to extract this information. 

5 Conclusions and outlook 

We have compared NLO predictions with data of the E706 experiment for various observables. 
We have classified the observables according to their respective infrared sensitivity. For the class 
of infrared insensitive observables the NLO predictions agree rather well with the corresponding 
measured ones over the whole measured spectrum. In contrast, the NLO result is inadequate to 
describe infrared sensitive observables, not only in the close vicinity of the infrared sensitive region as 
expected, but even over a large fraction of, if not the whole spectrum. The phenomenological analysis 
performed in this article suggests that the discrepancy can be understood in terms of a probabilistic 
reshuffle of the partonic configurations involved in the NLO calculation. From a theoretical point of 
view a Sudakov-type summation of logarithmically enhanced infrared terms which appear order by 
order in perturbation theory in the expansion of these infrared sensitive observables is mandatory. 
Yet it is likely to be not enough, as the discrepancy is probably dominated by large non perturbative 
effects. The non perturbative dynamics strongly distorts the infrared sensitive observables because 
the experimental selection cut pTmin imposed on each pion is not large compared to the typical 
non perturbative scale 0(1) GeV. On the other hand, these non perturbative effects turn out to be 
weak on infrared insensitive distributions. In particular they are much weaker than in the case of 
one particle inclusive production, especially on the px distribution in the lower px range. This can 
be traced back to the absence of trigger bias in pair production when symmetric cuts are applied, 
contrarily to the situation in one particle production. A further examination of the E7G6 data tends 
to suggest that these non perturbative effects become less and less important in the infrared sensitive 
distributions, and in particular more and more confined to the vicinity of the critical points, as the 
Pt's of the selected pions increase for fixed \/S, although, to be fair, no definite conclusion can be 
extracted from this examination due to large experimental uncertainties. On the other hand no 
evolution of the magnitude of non perturbative effects can be seen when going from 530 to 800 GeV 
because the level arm in ^/S is not large enough. 

We expect these features to hold as we move from fixed target to collider experiments, as soon as 
the selection cut pTmin increases up to 10-15 GeV at the Fermilab Tevatron and 20-40 GeV at the 
forthcoming LHG. We therefore expect NLO predictions to be reliable in particular for inclusive 
observables concerning pions pairs - as well as 777° pairs - at the Tevatron and LHG, at least 
far enough from the boundaries of the spectrum. On the other hand, near these boundaries, the 
various kinematic cuts can induce extra infrared sensitivity in partonic calculations. This is certainly 
expected to hold for observables which do not require any isolation cut. On the other hand, the 
implementation of isolation can affect the ability of fixed order perturbative calculations to provide 
reliable predictions for most inclusive distributions. This problem, which is of particular relevance 
in the calculation of backgrounds to the Higgs boson search at LHG in the channel if ^ 77, will be 
examined in detail in a forthcoming article ||3^ . 
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Figure 1: Scale dependence of the NLO prediction for the pr distribution in 800 GeV proton collisions on 

H target, /jP = (pri +Pt2)/2. 
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Figure 2: Dipion differential cross section da/dpr vs. pr, the transverse energy of each pion, in p — Be 
collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [^. The NLO prediction with scales M = /x = Mf = 3/8 (pri + PT2) is 
shown as triangles. 
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Figure 3: Dipion differential cross section da/dm-mr vs. m^^, the invariant mass of the pair, in p — Be 
collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration The NLO prediction with scales M — /i — Mf =3/8 (pri + PT2) is 
shown as triangles. With the chosen ordinate scale in the bottom histogram, the first two points fall outside 
the plot. 
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Figure 4: Dipion difFerential cross section da/dyTTTr vs. j/^^, the rapidity of the pair in p — Be collisions with 
a beam energy of 530 GeV. Here da/dyTTTr stands for J dm.nw{da / dy-n-ndm-mr) integrated over the following 
invariant mass slices: 5.0 < m,r7r < 6.5 GeV (top left), 6.5 < m-^-^ < 6.5 GeV (top right), 8.0 < m^^ < 10.0 
GeV (bottom left) and 10.0 < m^^ < 14.0 GeV (bottom right). Data points with statistical and systematic 
errors in quadrature are from the E706 collaboration [Q]. The NLO prediction with scales M — fi — Mf — 
3/8 (pTi + Pt2) is shown as triangles. 
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Figure 5: Dipion differential cross section da/dcos 9* vs. the variable cos 9* defined in subsect. 3.1, in p—Be 
collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [0. The NLO prediction with scales M — = Mf = 3/8 (pri + PT2) is 
shown as triangles. 
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Figure 6: Dipion differential cross section dcr/d0,rir vs. 0,r7r, the azimutlial angle difference between the 
two pions in a pair, in p — Be collisions with a beam energy of 530 GeV. Data points with statistical 
and systematic errors in quadrature are from the E706 collaboration The NLO prediction with scales 
M = /i = Ml =3/8 ijpTi + Pt2) is shown as triangles. 
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Figure 7: Dipion differential cross section dcr/dpT.out vs. the variable pT.out defined in subsect. 3.2 



m p- 



-Be 



collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [0. The NLO prediction with scales M — = Mf = 3/8 (pri + PT2) is 
shown as triangles. 
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Figure 8: Dipion differential cross section da/dqr vs. qx, the transverse momentum of a pair, in p — Be 
collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [0. The NLO prediction with scales M — = Mf = 3/8 (pri + PT2) is 
shown as triangles. 
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Figure 9: Dipion differential cross section da/dZ vs. the variable Z defined in subsect. 3.i, in p — Be 
collisions with a beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [0. The NLO prediction with scales M — fj, — Mf =3/8 {pri + PT2) is 
shown as triangles. 
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Figure 10: Ratios Rf{z) = D^/Y{z,Mf)/D^^j>'^ {z,Mf) of BKK to KKP fragmentations functions, vs. 
for various partons flavours /, at Mf = 3 GeV. Top: f = u quark [f = d quark equal to u). Bottom: / 
gluon. 
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Figure 11: Collisions of 530 GeV protons on Be target. Comparisons of E706 data with NLO predictions 
for the one pion inclusive pt spectrum, for BKK and KKP sets of fragmentation functions. The scale choice 
used is /i = M = M/ = pt/3. Only data points with pr > 4.35 GcV arc kept in order that the fragmentation 
scale Mf in the NLO predictions is above the starting scale %/2 GeV of the BKK and KKP sets. 
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Figure 13: Illustration of trigger bias effect in one particle inclusive production. The large circle has radius 
PTmin- The thick arrow stands for the transverse momentum of the pion generated by the hard subprocess. 
Left: pion not hard enough and no kr kick: event discarded . The radius of the inner circle on the right 
indicates the magnitude of the kr kick. Right: the pion would not be produced hard enough, but the 
partonic subprocess is kicked transversally: event accepted. 
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Figure 14: Illustration of the absence of trigger bias effect in inclusive pair production. The large circle 
has radius pr min ■ The thick arrows stand for the transverse momenta of the pions generated by the hard 

subproccss. The radius of the inner circle on the right indicates the magnitude of the kr kick. Top: the 
pions arc not produced hard enough by the partonic subproccss. Top left: no kx kick: event discarded. Top 
right: one pion is favoured by fcr kick, the other one is penalized: event also discarded. Bottom: the pions 
would be produced hard enough by the subproccss. Bottom left: no kr kick, event accepted. Bottom right: 
one pion is kicked down, even discarded. 
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Figure 15: Illustration of the reshuffle of back to back configurations by final state kr kick. The radii of 
the circles represent the magnitude of the kick. If the magnitude of the kr kick is not much smaller than 
the selection cut pr min , the 0,r7r of the reshuffled configuration can be very different from n. 



„ 400 

01 
CD 



A 

K 

-e- 
V 



300 



150 - 



100 



O E706 Be 530 GeV 
▼ NIC theory 



4^ 



10 11 12 13 14 

m„„ [ GeV ] 



Figure 16: Dipion distribution < ci(/>,r7r > /dm^.^ vs. the invariant mass m^^, in p — Be collisions with a 
beam energy of 530 GeV. Data points with statistical and systematic errors in quadrature are from the E706 
collaboration [Q]. The NLO prediction with scales M = /i = Mf =3/8 (pti + Pti) is shown as triangles. 
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Figure 17: Dipion differential cross section Dipion differential cross section da / dpr out vs. pTout in p — Be 
collisions with a beam energy of 530 GeV. Here da/dpTout stands for J dmTr-w{da / dpr outdrriTyTy) integrated 
over the following invariant mass slices: 5.0 < m^^ < 6.5 GeV (top left), 6.5 < mTrvr < 6.5 GeV (top 
right), 8.0 < m,r7r < 10.0 GeV (bottom left) and 10.0 < m^^ < 14.0 GeV (bottom right). Data points with 
statistical and systematic errors in quadrature are from the E706 collaboration ||^. The NLO prediction 
with scales M = jj, = Mf = 3/8 (pri + pt2) is shown as triangles. 
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Figure 18: Dipion differential cross section da/dZ vs. Z in p — Be collisions with a beam energy of 
530 GeV. Here da/dZ stands for J drnTrTrida / dZdniTrTr) integrated over the following invariant mass slices: 
5.0 < Tn-jrTT < 6.5 GeV (top left), 6.5 < m^^ < 6.5 GeV (top right), 8.0 < m^-,r < 10.0 GeV (bottom left) and 
10.0 < rriTrTr < 14.0 GeV (bottom right). Data points with statistical and systematic errors in quadrature 
are from the E706 collaboration [0. The NLO prediction with scales M = fj, — Mj — 3/8 (pti + Pt'/) is 
shown as triangles. 
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